Electronic Correlation and Transport Properties of Nuclear Fuel Materials 
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The electronic structures and transport properties of a series of actinide mono carbides, mono 
nitrides and dioxides are studied systematically using a combination of density functional theory and 
dynamical mean field theory. The studied materials present different electronic correlation strength 
and degree of localization of 5 /-electrons, where a metal-insulator boundary naturally lies within. 
We also investigate the interplay between electron-electron and electron-phonon interactions, both 
of which are responsible for the transport in the metallic compounds. Our results allow us to gain 
insights into the roles played by different scattering mechanisms, and suggest how to improve their 
thermal conductivities. 



When engineering fuel materials for nuclear power, im- 
portant thermophysical properties to be considered are 
melting point and thermal conductivity. Uranium and 
plutonium oxide fuels used in very high temperature fast 
breeder reactors have very high melting points, but they 
suffer from poor thermal conductivity, because in these 
insulating oxides only lattice vibrations conduct heat. 
Hence attention is turning to metallic fuels for the new 
generation of reactors, such as uranium carbide and ni- 
tride [H. Understanding the physics underlying trans- 
port phenomena due to electrons and lattice vibrations 
in actinide systems is a crucial step toward the design 
of better fuels. In this Letter we conduct a systematic 
theoretical study on the electronic structures and lattice 
dynamics of actinide compounds, and show that systems 
close to the Mott transition from the metallic side are 
the best option. 

Thermophysical properties of solids are determined 
from their electronic structures, but in actinides they are 
not well described by the traditional approaches based 
on density functional theory (DFT) within its local den- 
sity approximation (LDA) due to strong electronic cor- 
relation. In this study we use an advanced electronic 
structure method based on the combination of LDA and 
dynamical mean field theory (LDA+DMFT) Q, which 
has proven success in describing such strongly correlated 
problems 0-Q.Our full-potential charge self-consistent 
implementation of LDA+DMFT described in Ref.[8| is 
based on the DFT program of Ref. [3] ■ For the impurity 
solver we use the continuous time quantum Monte Carlo 
(CTQMC) algorithm d OH- For late actinides such as 
Pu and beyond, we use the less expensive vertex corrected 
one-crossing approximation (OCA) [2], which is very ac- 
curate in these more localized systems. All calculations 
were performed in the paramagnetic phase. 

We first describe the chemical trends governing the 
degree of localization of the /-electrons in the binary 
actinide compounds listed in Fig. [TJ The key parame- 
ters are: the onsite Coulomb repulsion among the 5/- 
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FIG. 1. Correlation diagram. The shading represents the 
electronic correlation strength. The labels on the top denote 
the actinides elements, and the labels to the left denote the 
ligand elements. The red line is the metal-insulator boundary. 
Two quantities are listed in each cell: Hubbard U (units: eV) 
and /-electron valence n,f. 



electrons, quantified by the Hubbard U and Hund's rule 
exchange J; the charge transfer energy A = s f — s p ; and 
the 5/ band width quantified by the hybridization be- 
tween 5/ and spd electrons. The charge transfer energy 
increases vertically from carbides to oxides due to the 
change in the electro-negativity of ligand atoms. The 
band width of 5/-electrons shrinks horizontally from U 
to Cm compounds, indicating a more localized nature in 
late actinides. This causes a reduction of screening which 
is manifest in the gradual increase of U from the left to 
the right, and from the top to the bottom of the table. 

While most electronic structure methods can accu- 
rately calculate the hopping integrals between various 
electronic orbitals, evaluating the screened U in solids 
is generally a difficult task. Here we have computed U 
using a newly developed fully self-consistent many-body 
GW approach [ll|, which provides a seamless interface 
with LDA+DMFT. The latter method allows to deter- 
mine the degree of localization of the 5/-electrons in 
each material. Our estimates for the Hund's J are within 
the range of 0.5 — 0.6 eV, about 30% smaller than their 
atomic values due electronic screening. As a combination 
of the above quantities, the overall correlation strength 
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FIG. 2. Theoretical DOS compared with available x-ray pho- 
toemission spectroscopy, (a) Partial 5/ DOS of UC and select 
actinide nitrides. The XPS & BIS data of UC is from Ref. Ill, 
UN from Ref.Ql], PuN from Ref.QjJ, and AmN from Ref. [Hi] . 
(b)Total and partial DOS of U0 2 . XPS and BIS taken from 
(c) Total and partial DOS of Pu0 2 . XPS from 
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and localization is visualized by the shading of Fig. [T] re- 
ferred as the "correlation diagram" of binary nuclear fuel 
materials, where the gray gradient approximately repre- 
sents the partial / density of states at the Fermi level 
computed by LDA+DMFT. 

Next, we present the frequency dependence of the elec- 
tronic spectral functions of some representative com- 
pounds in Fig. Ufa). From the top panel to the bot- 
tom, the 5/ partial DOS changes qualitatively. UC and 
UN represent an itinerant 5/-electron system with most 
spectral weight on the Fermi level, but the picture starts 
to change at PuN, where the Kondo resonance and satel- 
lite 5/ states are present. In AmN the 5/ DOS begins 
to form an marginal energy gap. The evolution of the 
density of states from UN to CmN echoes the itinerancy- 
localization transition of 5 /-electrons, and demonstrates 
the metal-insulator transition in a transparent point of 
view. CmC, CmN, and all the actinide oxides are also 
found to be insulators. This allows to establish a metal- 
insulator transition boundary, illustrated by the red line 
in Fig. [U 

The actinides ions in most of the metallic crystals are 
found to be in a mixed valence state, where they do not 
settle in one valence, but fluctuate between different va- 
lences in the solid. It can be described by an effective 
number nt (listed in Fig. [T]), obtained using a valence 
histogram technique [f| , which represents an average over 
all the atomic configurations weighted by corresponding 



probabilities. 

The total and partial DOS of UO2 and Pu02 are shown 
in Fig. &[b) and[2Jc). Most noticeably, the situation 
U > A allows us to describe the insulating actinide ox- 
ides as charge transfer Mott-insulators 12j , which is well 
known from late transition metal oxides, for example 
NiO, the classical textbook example of strongly corre- 
lated systems. Thus they could exhibit the Zhang-Rice 
state (ZRS) [l3j], which is the low-energy resonance cor- 
responding to the coupling of local moments of correlated 
electron orbitals to the hole induced by phototemission 
process on ligand orbitals. The appearance of the ZRS 
in UO2 can be understood from the existence of a local 
magnetic moment in the U 5/ 2 many-body ground state, 
which is the triplet. It is also clear that Pu02 does not 
show the ZRS because its ground state of the 5/ 4 -shell 
is the Ti singlet. 

Experimentally, UC is a Fermi liquid (FL) at room 
temperature and ARPES measurement indicates that the 
overall band width is reduced by a factor of 4 relative to 
the LDA band structure [20| . In our calculation UC is a 
FL below 300K with m*/m LDA = 3.7. On the other 
hand, UN shows a strongly correlated heavy fermion 
character with a coherence temperature below its Neel 
temperature of 53K. In the absence of magnetic order, 
UN would be a FL at very low temperature with a large 
mass enhancement (m* /jtilda ~ 12) as can be inferred 
from the linear specific heat coefficient It is a non 
FL in the temperature range (55 — lOOOif ) we studied. 

After understanding the electronic structures, we turn 
to the transport properties. We focus on correlated 
metallic compounds, where electrons play the role of 
charge and heat transporters, while retaining a high melt- 
ing point. Although in normal metals electron-phonon 
scattering is dominant except at very low temperatures, 
in strongly correlated metals electron-electron scattering 
takes the lead. The electronic contribution to thermal 
conductivity is proportional to the electrical conductiv- 
ity via the Wiedemann-Franz law. From the electronic 
structure and correlation strength of the studied materi- 
als, small resistivity occurs in the least correlated com- 
pounds in our table. Indeed UC and UN are the best fuel 
materials in terms of their outstanding transport prop- 
erties. 

Strong Coulomb interactions among electrons can sub- 
stantially reduce the interaction between electrons and 
lattice vibrations [2l|. Hence the electron-phonon in- 
teraction (EPI) is usually weaker in strongly correlated 
materials, which might lead to smaller resistivity due to 
EPI. On the other hand, increasing electronic correla- 
tions leads to an increase in resistivity due to electron- 
electron scattering. Therefore neither extremely weak 
nor strong correlations is good from the perspective of 
minimizing resistivity. Deciding the optimal degree of 
correlation for the purpose of maximizing conductivity 
thus requires first -principle calculations. 
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FIG. 3. Phonon dispersions of UN plotted along three high- 
symmetry lines. Solid curves: Calculated phonon dispersion. 
Open circles: Experimental phonon excitations from Ref . [25l] . 



For the calculation of phonon spectra we used the den- 
sity functional linear response approach implemented in 
the LMTO basis [H,|23j]. The spin-orbit coupling effect 
is included fully relativistically. The phonon dispersion 
of UN along 3 high-symmetry directions is plotted in Fig. 
EH together with experimental data measured by neutron 
scattering 25]. UC carries very similar phonon disper- 
sions but slightly lower phonon energies. Despite appar- 
ent presence of correlation effects, excellent agreement 
is achieved with the local density approximation (LDA). 
Similar success of LDA in studying lattice dynamics of 
strongly correlated metallic systems have been reported 
earlier, for example in cuprates [26| , and recently iron 
pnictides (27} . 

Calculations of electron-phonon interactions and 



transport properties require, on the other hand, quasipar- 
ticle description of the one-electron spectra when evalu- 
ating Eliashberg and transport spectral functions by in- 
tegrating over the Fermi surfaces [24j ]. As a result, due to 
large mass enhancement, the straightforward LDA pro- 
cedure can produce wrong electron-phonon resistivities, 
which was indeed found in our calculation for UC where 
p{T)epi was overestimated by a factor of 3 compared to 
experiment. 

In order to evaluate the electron-phonon scattering in 
the presence of correlations, we develop a method that 
accounts for the effects from quasiparticle mass renor- 
malizations and spectral weight transfer by utilizing in- 
teracting Green's functions. The use of the pole interpo- 
lation of self-energy allows us to replace the non-linear 
(over energy) Dyson equation by a linear Schroedinger- 



like equation in extended subset of "pole states" [31 1. 
The advantage of this method is that the well developed 
machinery of standard electronic structure methods can 
be simply generalized to account for the dynamical self- 
energy effects. 

In the present work we extend this method to compute 
the electron-phonon interactions for UC and UN, whose 
/-electrons show itinerant behavior with m* jm^uA ~ 



4 — 12. To capture this mass renormalization effect, 
we first make a fit to the self-energy obtained from the 
CTQMC, using a two-pole interpolation where the slope 
of the self-energy at zero frequency dY>(uj) / duj^^o = 
1 — m* /rriLDA controls the electronic mass enhancement, 
and the positions of the two poles determine the transfer 
of the spectral weight from the quasiparticle band to the 
Hubard bands. Next, we assume that the /-electrons 
are rigidly bound to their ions so that there is no actual 
change in the self-energy, caused by ion displacements 
from their equilibrium positions. Since the main contri- 
bution to electronic transport comes from the states near 
the Fermi surface, where quasiparticles are best described 
in terms of slave bosons, the neglecting of 5T,(u) due to 
ion displacements corresponds to a rigid self-energy ap- 
proximation. This is very similar to the famous rigid 
muffin-tin approximation (RMTA) [28j], which has been 
successfully applied in the past to study electron-phonon 
interactions in transition metal compounds [HHoi. The 
electron-phonon scattering matrix element gkj,k+qj' is 

(e) 

then evaluated using the electronic components ipU that 
appear as solutions to the extended Schroedinger-like 
equation. These matrix elements can be subsequently 
used for evaluating the EPI part of electrical and thermal 
resistivity [24j . where the corresponding Fermi surface 
integrals are now performed with the extended "band 
structures" that acquire renormalizations due to cor- 
relations. Using this method, a reduction (by a factor 
of 3) in p(T) epi for UC was obtained, while in UN the 
effect was marginal. 

To evaluate the conductivity due to electron-electron 
scattering we use the Kubo formalism 8j , where the scat- 
tering rate comes from the imaginary part of DMFT self- 
energy £(w,T), obtained from CTQMC. 

We can now build the entire picture of the electronic 
transport in the uranium compounds with our results 
summarized in Fig. 2J Electron-electron scattering can 
account for approximately 80% of p(T) in UN, commonly 
found in heavy fermion systems, entitling it as a strongly 
correlated bad metal. In contrast, UC shows nearly linear 
p(T), which is an indication of dominant electron-phonon 
scattering, and our calculated results indeed show that in 
UC, p(T) ee is much smaller than p(T)epi- Our findings 
verify the distinct characters in the electrical transport of 
UC and UN, two seemingly similar materials. By making 
a side-to-side comparison, we see that UC is less corre- 
lated than UN, which makes p(T) ee larger in the latter. 
Nevertheless p(T)epi is larger in the former. 

We also estimate lattice vibrational contribution to 
thermal conductivity in UC and UN. This is done by 
evaluating the Gruneisen parameter and phonon group 
velocities using the method described previously for ox- 
ide fuels @. At T = lOOOif, lattice thermal conductiv- 
ity n p h = 2.7W/mK in UC, and n ph = AAW/mK in 
UN, much smaller than electronic counterpart. Thus K p h 
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FIG. 4. Electrical resistivity due to two different scattering 
mechanisms. The electron-phonon interaction (EPI) resistiv- 
ity is show as solid lines, and electron-electron interaction 
resistivity, which is calculated by LDA+DMFT (CTQMC) at 
several temperatures, is shown as stars, (a) UN. Experimental 
resistivity data are taken from Ref. [l5|]. (b) UC. Experimental 
data after Ref. [33-l3a] . 

only plays a minor role in these two metallic uranium 
compounds. 

We put together our results and evaluate total ther- 
mal conductivity at lOOOif , a representative tempera- 
ture under which nuclear reactors operate. By applying 
the Wiedemann-Franz law on the electrical conductivity 
data, we obtain n ee . Since electronic thermal resistivity 
consists of two scattering processes, total thermal con- 
ductivity is estimated by n to tai = (k^ 1 + ^p/)^ 1 + K ph, 
in which the first two terms correspond to ^electron- For 
UN, our result is n to tai — W.5W/mK, compares well 
with a recent study which extracted the phonon con- 
tribution from molecular dynamics (MD)[38jj and the 
electronic contribution from experiments. Experimen- 
tally, K(lOOOiT) w 19 - 23W/mK. In UC, we obtained 
ntotai = 18.7WJmK, also close to the experimental value 
of 23W/mK [33]. The discrepancy between theory and 
experiment is likely due to other excitations that can con- 
duct heat and are not accounted for in our calculation, 
and the approximate nature of the Wiedemann Franz 
law and Boltzmann transport theory which are used to 



obtain the electronic and lattice thermal conductivity, 
respectively. 

At last, the understanding gained from our compu- 
tational study suggests avenues for improving the ther- 
mal conductivity of UC and UN. Since optimizing ther- 
mal conductivity is equivalent to minimizing resistivity at 
high temperatures, it is interesting to look at the doping 
dependence of the resistivity in UC/UN, or in the solid 
solution UCi-zN^. Let us represent the total resistivity 
of UCi-.N, by p(T,x) = p(T,x) ee + p(T,x) E pi- Note 
that on one hand, p(T,x) ee should show a rapid growth 
when x approaches 1 since p(T) ee ex 1/Tk oc W/Z above 
coherence temperature, where Tr is the Kondo temper- 
ature and W is the band width. On the other hand, 
p(T, x)epi being proportional to transport coupling con- 
stant Aj r divided by the average square of the Fermi ve- 
locity N(0)(v 2 ) of the quasiparticles decreases with x in 
our calculation. Therefore there exists a region between 
UC and UN where the total resistivity is minimized. It is 
also possible to achieve similar effects in UC by electron 
doping, or in UN by hole doping. 

To conclude, we have carried out the first 
LDA+DMFT exploration of the electronic struc- 
ture and transport properties of nuclear fuels. The 
actinidc dioxides are charge-transfer insulators, where 
the Zhang-Rice state is present in UO2. The metallic 
carbide and nitride compounds exhibit strong electronic 
correlations, which is reflected in the incoherent non 
Fermi liquid behavior at temperatures relevant for 
nuclear reactions. We have achieved a successful theo- 
retical description of the transport in UC and UN. While 
UN clearly shows a strongly correlated signature, both 
the electron-electron and electron-phonon scattering 
mechanisms contribute to transport in the less correlated 
sister compound UC. Our findings enable us to give 
predictions on how to improve these two uranium based 
nuclear fuel materials. 
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